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We consider the effects of certain forms of decolierence applied to botii adiabatic and non-adiabatic 
geometric phase quantum gates. For a single qubit we illustrate path-dependent sensitivity to 
anisotropic noise and for two qubits we quantify the loss of entanglement as a function of decoherence. 
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The discoveries of quantum algorithms for factoriza- 
tion [0 and searching , and techniques for quantum er- 
ror correction j^, ^ (and subsequent fault-tolerant meth- 
ods) have generated considerable motivation for the re- 
alisation of quantum computing (QC) hardware. Signif- 
icant progress has beerr made at the few qubit level; at 
present many possible alternative routes are under ex- 
ploration Q. The use of fundamental entities as qubits 
currently leads the way; however, there is growing longer 
term hope that fabricated condensed matter systems 
(maybe still using fundamental rather than fabricated 
qubits) may provide the vehicle for scalability in qubit 
number. Nevertheless, at present the decoherence prob- 
lems in such systems loom very large indeed. Fault- 
tolerant techniques cannot be brought into play unless 
the underlying decoherence rates are small in the first 
place. 

As arbitrary single-qubit gates and some entangling 
two-qubit gate are universal for quantum computing 
1^, 0, gl, experimental QC focuses on reahsing such gates. 
Whilst it would be unfair to deem single-qubit gates triv- 
ial, the entangling of qubits represents the first major 
hurdle — one cannot claim to have a serious QC candi- 
date until this has been achieved. Nevertheless, with any 
new contender it is natural to investigate the simplest 
gates first, before progressing to those based on qubit 
coupling. Clearly a method for inferring the likely level 
of entanglement in a two-qubit gate from the results of 
single qubit experiments is a handy tool. This is the basis 
of the simulation approach we present here in the form 
of a specific example. It is certainly true that detailed 
calculations of decoherence effects (based on tracing over 
microscopic environments) can yield considerable under- 
standing of the form of decoherence seen by a qubit, and 
can sometimes give estimates of decoherence rates; how- 
ever, in any experimental realisation the ultimate calibra- 
tion of decoherence is through measurement. Given this, 
we consider a simple simulation approach which, based 
on the observations of single qubit gates, enables estima- 



tion of the level of entanglement (if any!) that could be 
expected in a two-qubit gate with similar qubits. This re- 
lies on some advanced knowledge of the form of the dom- 
inant decohering mechanisms (although some of this may 
be inferred from various independent single-qubit exper- 
iments) and an assumption of environment effects in in- 
teracting two-qubit experiments being similar to those in 
the single-qubit experiments. If the former are worse, for 
example due to additional external source terms in the 
Hamiltonian, then an optimistic upper bound on entan- 
glement results. Despite these constraints, the simulation 
approach could be very useful for new experimental QC 
investigations. 

To illustrate the approach, we apply it to geometric 
phase gates H, ^ 0. This provides an example of the 
approach; however, the specific case of geometric phase 
gates is of interest in its own right, since the technique 
has already been applied to NMR experiments [^, has 
been proposed for use with superconducting qubits 
and, in principle, can be applied to other realisations of 
qubits with suitable source terms in their Hamiltonians. 
We quantify the effects of different forms of decoherence 
applied to single qubit phase gates, and the loss of en- 
tanglement for a conditional two-qubit gate. The same 
approach can also be applied to all forms of dynamically 
generated gates; a broader detailed study will be pre- 
sented in a future paper. 

The model of decoherence we use is Markovian, with 
the reduced density operator p of the qubit system de- 
scribed by a Bloch-type master equation 
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(1) 

where H is the system Hamiltonian, the operators {Lm} 
represent the coupling to the environment and h — I. 
The implicit origin of the non-unitary evolution gener- 
ated by the {Lm} is coupling to a bath of environment de- 
grees of freedom, which are traced out to give the reduced 
p. The form of Eq.(|l|) is somewhat restrictive, but within 
this Markovian limit it is possible to describe phenomena 
such as dissipation (spontaneous decay and finite temper- 
ature stimulated effects), white noise Hamiltonian terms 
and quantum measurement interactions. It is therefore 



2 



possible to treat a number of realistic forms of decoher- 
ence. The master equation (|^) can be solved in many 
simple cases; however, in order to be able to treat the 
types of time dependent Hamiltonians (including pulses 
etc.) used for the realisation of actual quantum gates, we 
generally work numerically. We use a quantum trajectory 
method, quantum state diffusion jl^, |l3, 14 1, to solve the 
master equation (|l]) through averages over stochastically 
evolving quantum states. Whilst we only give statisti- 
cal data here, it is worth noting that (ensemble NMR 
systems aside) since actual quantum gates/computations 
run on individual systems, such simulation techniques in 
fact produce results in a manner very akin to actual QC 
experiments and examination of individual trajectories 
can provide additional insight |jl^. In all examples pre- 
sented here the averages are over 1000 trajectories unless 
otherwise stated. 

Our first example is a single qubit geometric phase 
gate. In order to generate a purely geometric phase, the 
dynamical phases acquired by the different amplitudes in 
an arbitrary qubit state have to be cancelled ||l^, . We 
use the scenario of Ref.lO. The path traversed in param- 
eter space is amenable to the investigation of different 
forms of decoherence and it has also been implemented 
experimentally [^. A spin qubit (| tz) = |0) , | iz) = |1)) 
is subject to a static z-magnetic field ujo and (within the 
usual rotating wave approximation) a field of amplitude 
LUi and phase 4> at angular frequency lu. The general 
Hamiltonian is 

H — -^<^z + [cos{Lot + <p)<7x + sin{uot + 0)cry] , (2) 

and to effect an ideal phase gate the spin is subjected to 
the unitary sequence U^^ = Iif C TUT C T . Here T is 
a tipping of the magnetic field through angle 6 (ramping 
uji from zero, with cosO — (wq ~ cj)/^ (wq — w)^ + lo\ 
with (j) at zero, C is a 27r rotation of the phase (j) 
fixed LUi and the bars denote the reversed paths. These 
operations have to be carried out adiabatically to avoid 
errors in the spin component amplitudes [ pO[ . Fast tt- 
pulses n interchange the t smd J, amplitudes half way 
through and at the end (to cancel the dynamical phase 
contributions). The ideal gate to effect a relative phase 
of 75 =47 on IV') = 2-i/2(|0) |1)) is 
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(exp(-2z7)|0) + exp(2i7)|l)), (3) 



where 7 = 7r(l — cosO) is the solid angle subtended by C 
at the origin. Part of the appeal of this form of quantum 
gate is its potentially different sensitivities to different 
forms of decoherence; a potential drawback is the need 
for adiabaticity, so the gate is slow and decoherence has 
longer to bite. 

We have studied these effects in detail. Results are 
shown in Fig. 1(a) for the effects of noise in the x- or 
z-components of the magnetic field, for two different 75. 
The gate was run adiabatically giving a zero decoherence 



fidelity Q of / = 0.999984 for 7^ = tt and / = 0.999993 
for 7b = 7r/8. For the smaller 7^, which corresponds 
to a tipping of 14.3615 deg (so the instantaneous energy 
eigenstates remain closer to t and | in z), the system is 
clearly significantly more sensitive to z-noise compared 
to x-noise On the other hand, for the larger 73 

where the tipping is 41.4096 deg there is less distinction. 
The case of isotropic noise is illustrated in Fig. 1(b). 




0.005 0.01 0.015 0.02 0025 003 0.035 0.0+ 



5 0.02 0.025 003 0.035 004 

r 



FIG. 1: The effects of anisotropic noise (generated by L = hia^ 
or L = KfJz) on the gate Eq.dS) are shown in a) for loq — lo = 
100. For adiabaticity T takes (dimensionless) time tt and 
C takes 2tt. The II's are square pulses taking time tt/IOO. 
The ratio of these slow and fast operations is comparable 
to the frequency ratio in The loss in fidelity (1 — /) is 
shown against the noise parameter F = k^: 7s — tv/8, z- 
noise (i) and a;-noise (iv); 7s = tt, z-noise (ii) and a;-noise 
(iii). In b) the effects of isotropic noise (generated by three 
independent Li = Kai for i — x,y,z) on the gate Eq.(|^ for 
7s = TT. (1 — /) (i) and entropy S = —Tr{p\og2 p) (ii) are 
shown against F = k^. 

The small-F rate of fidelity loss is twice that of the worst 
behaviour in Fig. 1(a) (which follows from a simple ana- 
lytic estimate) and indeed the whole fidelity loss fits well 
with the analytic form 1 — / = i(l — exp(— 4Fr)) where 
T is the gate duration. Provided that the gate is adia- 
batic, the effects of isotropic noise are set by the gate 
length and level of decoherence, independent of the gate 
details. Also shown in Fig. 1(b) is the final system en- 
tropy as a function of F. For small F the entropy (loss 
of information) increases significantly faster than the loss 
in fidelity 

Our second example is the more important case of a 
conditional two-qubit geometric gate [|lOj, where entan- 
glement is generated, or not, depending upon the level 
of decoherence. This requires two spin qubits with bare 
transition frequencies of uja (target qubit) and Wf, (con- 
trol qubit) and coa 3> uJb- An interaction Hamiltonian of 
Hmt = jcraz<^bz generates the conditional phase. This 
form of interaction is that appropriate for NMR and 
certain condensed matter qubits. We have left J fixed 
for the simulations presented here (as is appropriate for 
NMR systems), but in principle this coupling may be 
tunable for some condensed matter scenarios. A condi- 
tional phase gate (introducing A7 = 7_b(T6) — 7B(ib), for 
the two states of the control) is realised by the unitary 
sequence Ua-j = H^C/^^H^C/"^ . Here superscripts refer 
to the qubit operated upon and U^^ is the same as U^j 
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but with the final 11 removed. For \'i')ab = |'0)aIV')6 this 
generates 

UA^mab = \ (cxp(-2iA7)|0),|0)b + cxp(2^A7)|0)„|l)^, 

+ exp(2iA7)|l)a|0)6 + exp(-2iA7)|l)a|l)6) .(4) 

For a phase of A7 — j this state has a concurrence 
of unity and so is a maximally entangled two-qubit state. 

We have chosen the coupling and frequency parameters 
so that the zero decoherence state at the end of the simu- 
lated conditional phase gate is maximally entangled with 
fidelity of 0.999946, and then investigated this gate un- 
der various forms of decoherence. Examples of the results 
are shown in Fig. 2. Clearly the rate of fidelity loss and 
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FIG. 2: The effects of equal isotropic noise (Li = kct^ for 
i = x,y,z) applied to both qubits for the gate of Eq.(^ are 
displayed in a) with A7 — ^. Shown are plots of the loss 
in fidelity (i), the entropy S = — rr(plog4p) (ii) and the 
entanglement of formation (EOF) (iii) as a function of F = . 
The gate timings are as in Fig. 1 and additional parameters 
used are cOa — u) = 100, ujb = 1, Wi = 87.9238 and J = 37.5. 
The entropy is base 4 to ensure a maximum value of 1. Also 
shown in a) is an enlargement of the region near F = 0. In 
b) we also show the decrease in EOF as a function of F for Sx 
(i), Sz (ii) and isotropic (iii) noise on each qubit. 

the sympathetic increase in entropy are correspondingly 
greater for this system as noise is acting independently 
on the two qubits. For this system p was reconstructed 
by tomography [|l6| from the sixteen expectation values 
(o'aiO'bj) for i,j = 0, a;, y, z (ctq = Identity), akin to what 
is needed in any two-qubit experiment for a full recon- 
struction of p. From this it is possible to compute the 
entropy and some measure of entanglement. For illus- 
tration we use the entanglement of formation (EOF) 
as this gives an upper bound on the level of decoher- 
ence for which some entanglement can be said to exist 
i^thres = 0.00445 for isotropic noise in our example). 

The maximally entangling parameters used for Fig. 2 
generate relatively large tipping angles, so there is only a 
slight sensitivity to the direction of anisotropic noise ap- 
plied to both qubits, as shown in Fig. 2(b). Furthermore, 
detailed studies not illustrated here show that there is 
only a minor difference between the separate effects of 
noise on the control and the target, so from an exper- 
imental perspective on such a gate, there is nothing to 
be gained by singling out either of these for decoherence 



reduction measures (e.g. error correction). Both these 
points hold right down into the very small decoherence 
regime, where any practical system would have to oper- 
ate. 

In our simulations so far the gate times have been set 
to ensure adiabaticity. Such gates are relatively slow 
and therefore exposed to the ravages of decoherence for 
longer. Conventional dynamic gates can run much more 
quickly and so for comparison we investigate a dynamic 
gate based on the same interaction as the entangling adi- 
abatic gate. The gate is realised by the unitary evolution 
U{T) = exp \i^(Jaz<^bz\ where we again choose J = 37.5 
and T is now the interaction time. In the absence of de- 
coherence this produces the maximally entangled state 

|*)ab = ^e-*f (|0)a|0)fa + ma\l)b + l\l)a\Q)b + \l)a\l)b) 

with a fidelity of 1 if the gate acts for a total time Ttot — 
27r/75. For the adiabatic geometric phase gate Ttot = 
12.00047r, so the dynamic gate is approximately 450 times 
faster and its speed is limited directly by the strength of 
J. The results for this dynamic gate are illustrated in 
Fig. 3(a). As isotropic noise acts on the dynamic gate. 





FIG. 3: The effects of equal isotropic noise applied to both 
qubits for the dynamic gate are displayed in a) and to the 
nonadiabatic geometric gate in b). Shown in the plots are 
the loss in fidelity (i), the entropy S — —Tr{plog^ p) (ii) 
and the entanglement of formation (iii) as a function of F. 
For the dynamic gate the following parameters were chosen: 
J = 37.5 and T = ir/J. For the nonadiabatic fast geometric 
gate J = 37.5, Slo = 18.75, lob = 0.01. 

the entanglement of formation falls to zero at Tthres ~ 2 
(compared with Tthres = 0.00445 for the adiabatic gate). 
In fact we find that Tthres Ttot ~ An/lb for both the 
adiabatic geometric gate and the dynamic gate. There- 
fore, as the adiabatic gate operates for a significantly 
longer time, it is much more severely affected by deco- 
herence. This has serious implications for the physical re- 
alisation of such a geometric gate. Recently, however, it 
has been proposed to use the non-adiabatic, or Aharonov- 
Anandan, phase to speed up geometric phase gates p8| . 
Here the achievable reduction in decoherence is not en- 
tirely clear as this technique also introduces new poten- 
tial decoherence sources [|l^. We implement the uni- 
tary sequence i?S(7r/2)i?J(-7r/2);7(7r/4)i?°(-7r/2) where 
C7(^/4) = R-{3TT/4)U{TT/J)R-{-n/2)U{n/J)R-{-n/4} 
to produce the maximally entangled state \'^)ab- Here 
R%{0) indicates a rotation of about the x axis. The re- 
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suits of isotropic noise acting on both qubits are displayed 
in Fig. 3(b). As this fast geometric gate has two periods 
of free evolution it acts for twice as long as the dynamic 
gate. Hence, it is subject to the effects of decoherence 
for longer and entanglement is again lost at a faster rate 
(sHghtly greater than twice the rate, Tthres ~ 0.945). 
Overall, our results suggest that geometric gates proba- 
bly only offer a real advantage if they can be implemented 
faster than the equivalent dynamical gates. At present, 
such proposed gates do not beat the dynamic gate time 
(tt/J for our example), and it is not obvious that this 
is possible given both approaches realise entanglement 
through interacting qubit evolution. However, whether 
non-adiabatic geometric gates can be implemented more 
quickly is an open and critically important question cur- 
rently under investigation. 

A number of comments can be made in conclusion. 
(1) Single qubit geometric phase gates can show some 
sensitivity to the direction of anisotropic noise, but this 
is path-dependent. In our example, if the path is cho- 
sen to generate a large (~ it) relative phase between the 
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qub |t amplitudes, there is very little sensitivity. 



(2) From the general perspective of our simulation ap- 
proach, single qubit gate behaviour can be used with ex- 
perimental results to calibrate the level of decoherence 
present in a system p3|. Although the anisotropic noise- 
sensitive gates may be of limited use for actual quantum 
information processing (due to the small phase difference 
generated), they may be applied, for example coupled 
with the ability to reorientate the external static mag- 
netic field (z-axis), in mapping out the forms of decoher- 
ence acting on a qubit, in addition to calibrating them. 
This could be extremely useful for new experimental sys- 
tems where the dominant environmental coupling is un- 
clear in advance. 

(3) The single qubit decoherence calibrations can be used 



to p redict the expected level of entanglement in two-qubit 



gates (as illustrated in Fig. 2) prior to experiment. 
(4) The adiabatic gates discussed here are slow (relative 
to the timescale for dynamic gates) and so exposed to 
the ravages of decoherence for longer. To overcome this 
it is necessary to perform geometric gates faster than the 
equivalent dynamic gate. Whether they can be made 
faster than the dynamic gate is a unanswered question 
left for future investigation. The practical use of such 
geometric gates depends upon the resulting answer. 
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In principle the central T and T can be omitted if 11 is 
done by a TT rotation about the i/-axis (and we do so), 
but these central tips were employed in the experiments 
reported in [^. 

In terms of measurable quantities the gate fidelity is / = 
i [1 — sin6'cos7_B(o-z) -I- cos6'cos7_b(o-i:) -I- sin7s(o-y)]. 
In effect, there is more gate sensitivity to random fluc- 
tuations in the energy eigenvalues (the eigenstates stay 
close to those of a^) than to random precessions about 
the X-axis. 

All the quantities used in our simulations are 
dimensionless — an experimentally measured quantity 
such as a transition frequency sets the scale for every- 
thing in any actual qubit realisation. 
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